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(57) Abstract 

A system for locating objects in space, such as medical instruments (16) within the body of a patient (P), based upon transmission of 
magnetic fields from coils (10) in a fixed frame of reference to sensors (20, 22, 24) on the object or vice-versa. The current supplied to the 
coils is adjusted to assure that the sensors receive fields within a preselected range of magnitudes regardless of the location of the object 
in space. This assures mat the sensor operates within its optimum range, and permits use of compact transmitters and sensors. 
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DggCRIPTION 

MAGNETIC LOCATION SYSTEM 
5 WITH ADAPTIVE FEEDBAC K CONTROL 

TECHNICAL FIELD 

The present invention relates to systems for 
determining the location and/ or orientation of objects 
10 in space by detecting magnetic fields. 

BACKGROUND ART 

Various systems have been proposed for 
detecting the position and/ or orientation of an object 
using magnetic or electromagnetic fields. These systems 

15 typically employ field transmitters, such as 
electromagnet coils, disposed at known locations in a - 
fixed reference frame and a sensor, such as a coil or 
other transducer mounted to the object to be located. 
Each transmitter projects a field varying in space in a 

20 fixed frame of reference. The pattern of variation in 
space for each transmitter is different than the pattern 
for each other transmitter. For example the 

transmitters may be identical to one another but 
disposed at different locations or in different 

IS orientations. The field patterns of the transmitters 
are thus displaced or rotated relative to one another 
and relative to the fixed frame of reference. The sensor 
on the object detects the parameters of the field 
prevailing at the location of the object as, for 

JO example, the magnitude and/or direction of the field at 
the object or the magnitudes of individual components of 
the field at the object in one or more preselected 
directions. The transmitters may be actuated in a 
predetermined sequence so that at any time only one 

*5 transmitter is active and therefore the field prevailing 
at the object is only the field contributed by one 
transmitter, plus a background field due to the Earth's 
magnetic field and other environmental sources. 
Alternatively, the transmitters can be driven at 

K> different frequencies so that components of the signal 
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instruments such as endoscopes and catheters based upon 
transmitted fields are disclosed in United States Patent 
No. 5,042,486; 5,099,845; 5,211,165; 5,251,635; 
5,253,647; 5,255,680; 5,265,610 and 5,391,199. 

Systems of this nature typically have used 
fields having strength which varies as the third power, 
or a higher power, of distance from the transmitter. 
There is, accordingly, a very large variation in the 
strength of each field from location to location. In a 
medical device locating system, adapted to detect the 
location of the sensor anywhere within a sensing volume 
having dimensions on the order of 0.5 to 1.0 m, the 
field strength may vary through many orders of magnitude 
from one end of the sensing volume to another. When the 
sensor and the object to be located happen to be near a 
particular transmitter, they will receive an extremely 
strong field, whereas when the sensor and object are at 
the end of the sensing volume remote from the particular 
transmitter, they will receive only a very weak field 
from that transmitter. This requires a sensor with a 
very substantial dynamic range, capable of accurately 
monitoring very strong and very weak fields. That, in 
turn, poses stringent demands on the design of the 
sensor and may require a larger sensor. Moreover, the 
signal to noise ratio is poor when monitoring a very 
weak field and the accuracy of the system is accordingly 
less than optimum in these conditions. Both of these 
considerations are particularly important in systems for 
medical applications, where the sensor must be small to 
fit inside a medical device such as a catheter or 
endoscope. Although such systems can provide useful 
results, further improvement would be desirable. The 
same problems arise in systems using a coil or other 
transmitter on the object to be located and plural 
sensors in the fixed frame of reference. Here again, a 
particular sensor will be exposed either to a very 
strong field or a very weak field depending on whether 
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the sensor is at an unknown location within the sensing 
volume, and to provide one or more sensor signals 
representing such detected parameters. For example, the 
sensor may include a sensor body and a plurality of 
5 component sensors disposed on the sensor body, each 
component sensor being operative to measure the 
magnitude of a magnetic field component in a preselected 
local direction relative to the sensor body to provide a 
component sensor signal representing that particular 

10 component of the field prevailing at the sensor. The 
apparatus further includes calculation means for 
calculating the location of the sensor, the orientation 
of the sensor or both based upon the sensor signals and 
upon the unknown patterns of variation of the fields. 

15 The apparatus also includes feedback control means for 
adjusting the field generating means to alter the known 
pattern of variation of at least one of the fields 
responsive to the sensor signals, to the calculated 
location of the sensor or both so as to maintain the 

20 detected parameters of the altered field at the sensor 
within a preselected range. For example, where the 
field generating means includes a plurality of 
transmitters disposed adjacent the sensing volume, the 
feedback control means can be arranged to increase the 

25 strength of the field transmitted by a particular 
transmitter when the sensor is remote from such 
transmitter, and to decrease the strength of the field 
emitted by that transmitter when the sensor is close to 
the transmitter. Thus, the parameters of the field 

30 detected at the location of the sensor will always lie 
within a relatively narrow range of values. Because the 
transmitted field is altered in a known manner,, the 
position and/or orientation of the sensor can still be 
calculated, simply by using the altered field as the 

35 basis for calculation. For example, where the feedback 
control means increases the current flowing to an 
electromagnetic coil, the parameters used in the 
calculation to represent the strength of the field 
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a sensing volume. The apparatus according to this 
aspect of the invention includes calculation means for 
determining a calculated location of the transmitter 
based upon signals from the sensors and upon the known 
5 pattern of variation of the field at the transmitter. 
Here again, feedback control means are provided for 
adjusting the field generating means to alter the known 
pattern of variation of field sent by the at least one 
transmitter in response to the sensor signal, to the 

10 calculated location of the transmitter or both so as to 
maintain the parameters of the field from the at least 
one transmitter prevailing at each sensor within a 
preselected range. This arrangement is similar to the 
arrangement discussed above, except that the transmitter 

15 or transmitters are disposed on the object to be located 
whereas the sensors are located in the fixed frame of 
reference of the sensing volume. Here again, when the 
transmitter is close to the sensor, or where the 
orientations of the transmitter and sensor favor strong 

20 coupling between the transmitted field and the sensor, 
the strength of the field is reduced. If the 
transmitter is relatively far from the sensor, or in a 
unfavorable orientation, the control means increases the 
strength of the transmitted field. This arrangement 

25 provides advantages similar to those discussed above. 

A further aspect of the invention provides 
apparatus for sensing the disposition of an object — its 
position, orientation or both — in a frame of reference. 
Apparatus according to this aspect of the invention 

30 includes transmitter means including at least one 
transmitter for providing a field and sensor means 
including at least one sensor for sensing one or more 
parameters of a field prevailing at each such sensor and 
providing one or more sensor signals indicative of such 

35 parameters. The transmitter means and said sensor means 
cooperatively define a plurality of transmitter-sensor 
pairs, each including one transmitter and one sensor as 
elements of the pair, one element of each such pair 
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perspective view depicting portions of apparatus in 
accordance with a further embodiment of the invention. 

BEST MODE OP CARRYING OUT INVENTION 

Apparatus in accordance with one embodiment of 
5 the present invention includes three generally helical 
transmitter coils 10 disposed in a common plane. Coils 
10 are mounted in fixed position in the frame of 
reference of patient receiving bed 12. That frame of 
reference is denoted by a Cartesian coordinate system X, 
10 Y, Z as shown in Fig. l. A patient P may be positioned 
on the patient receiving bed. The axes 11 of the coils 
are parallel to one another. A patient receiving bed 12 
extends just above the plane of coil 10. The apparatus 
further includes an object or probe 14. The probe is 
adapted to be inserted into a medical instrument such as 
a catheter 16 and positioned at a desired location in 
the catheter, such as at the distal tip of the catheter 
or at another location along the length of the catheter. 
Probe 14 has mounted thereon a sensor 18. Sensor 18 
includes three component sensors 20, 22 and 24 adapted 
to sense components of the magnetic field in mutually 
orthogonal local directions X* , Y« and 2'. That is, 
component sensor 20 is sensitive to magnetic fields 
directed in direction X\ but largely insensitive to 
25 fields in directions V and Z', whereas sensor 22 is 
sensitive only to fields in direction Y' and sensor 24 
is sensitive to fields in direction Z'. These sensors 
are adapted to provide separate sensor signals 
representing the separate components. Sensor 18 may be 
30 a solid state sensor of the type described in the 
aforementioned International Patent Publication 
WO95/09562, the disclosure of which is hereby 
incorporated by reference herein. As further described 
therein, each of the component sensors may include a 
35 generally planar magnetically sensitive film, such as a 
magneto-resistive film or a Hall effect sensing film. 
Each such film may be sensitive to fields directed 
normal to the plane of the film. Alternatively, sensor 
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Coramand unit 28 further includes a control 
unit 34. Control unit 34 is linked by output lines 36 , 
38 and 40 to three separate coil drivers 42, 44 and 46. 
Each coil driver is linked to a separate one of 
5 transmitter coils 10. Each coil driver is adapted to 
send a direct current through the associated transmit 
coil 10. Each coil driver is arranged to control the 
amplitude of such current, and to turn the current on or 
off, in response to control signals received from 

10 control unit 34. The control unit is arranged to signal 
the coil drivers to provide currents to their respective 
transmit coils in alternating sequence, so that coil 10a 
receives current while coils 10b and 10c are inactive; 
coil 10b receives current while coils 10a and 10c are 

15 inactive and coil 10c receives current while coils 10a 
and 10b are inactive. The control unit receives data 
from amplification and conversion section 30 and, as 
further discussed below, actuates the coil drivers to 
vary the amplitude of the current to each coil. The 

20 control unit may include conventional interface devices 
such as digital to analog converters or bus interface 
units so that the output of the control unit is 
compatible with the control input of each coil driver. 
Also, although the control unit is illustrated 

25 separately from the other logical units of command unit 
28, it should be appreciated that the control unit may 
share physical elements of the command unit and other 
elements. For example, where the command unit 

incorporates a general purpose computer, the processor 

30 of the computer may serve both as an element of the 
position computation unit and as an element of the 
control unit, executing functions appropriate to the 
different units at different times. 

In a method according to one aspect of the 

35 invention, an catheter 16 is advanced into the body of a 
patient P. The probe 14 with sensor 18 thereon is 
disposed at the tip of the catheter. The catheter tip 
is disposed at an unknown location somewhere above the 
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all of the field magnitudes are within a preselected 
range of magnitudes. This preselected range is chosen 
to lie within the optimum operating range of sensor 18. 
Thus, the minimum field magnitude is selected to lie 
5 well above the noise threshold of the system and above 
the minimum sensitivity level of the sensor, where the 
maximum field level is selected to lie well below the 
maximum limit of linearity of the sensor, and well below 
the maximum field which the sensor can tolerate without 

10 loss of accuracy. For a typical magneto-resistive 
sensor, which is most accurate and repeatable when used 
with fields of less than about 4 Gauss, the preselected 
range of field magnitudes may be from about 1.0 to about 
2.5 Gauss. With typical Hall effect sensors, which are 

15 most accurate when used with fields over about 30 Gauss, 
the preselected range will be above about 30 Gauss. If 
the three field vector magnitudes observed during 
actuation of the three units 10a, 10b and 10c all lie 
within the preselected range, than the system proceeds 

20 to calculate the position and orientation of sensor 18, 
and hence the position and orientation of probe 14 and 
the associated tip 16 of the catheter using conventional 
position finding algorithms. For example, the 

mathematical methods disclosed in United States Patent 

25 4,710,708 for finding positions using a multiple 
transmitting or receiving stations and a multi-axis 
sensor can be employed. The disclosure of said f 708 
patent is hereby incorporated by reference herein. 
Briefly, the magnitude of the fields in each of the 

30 local or sensor directions X 1 , Y 1 , Z 1 represented by 
each of the component sensor signals from * each of the 
component sensors 20, 22 and 24 is a function of the 
overall strength of the field from the coil (also 
referred to as the magnetic dipole moment of the coil) , 

35 the distance from the particular coil to the sensor and 
the sensor rotation angles, i.e., the angles between the 
local directions X 1 , Y> and Z 1 and the coil frame of 
reference directions X, Y and Z. When three component 
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large change can be made when the field magnitude is far 
outside of the range, whereas a smaller change can be 
employed when the field magnitude is close to the range 
or within the range. The correction process continues 
until all field strengths are within the preselected 
range, whereupon the system calculates position and 
orientation. After the system finds values of coil 
currents which produce field magnitudes within the 
range, subsequent actuation cycles use these currents. 
During operation, as the physician utilizes catheter 16, 
the position of the catheter tip, and hence the position 
of probe 14 and sensor 18 may change. Any such change 
may bring one or more of the field magnitudes outside of 
the preselected range, whereupon the system will 
readjust the currents to the coils once again. 

As the feedback control unit readjusts the 
currents to the coils, the altered values of current are 
translated into new values for field strengths from the 
individual coils which are used in the aforementioned 
position determining equations. In this manner, the 
system assures that the sensor is always exposed to a 
field having a magnitude within the preselected range 
wherever the sensor is placed within a sensing volume 50 
extending over a preselected region above the plane of 
coils 10. The exact size of sensing volume 50 will 
depend upon the breadth of the preselected field 
magnitude range and the dynamic range of coil drivers 
42, 44 and 46, i.e., the degree to which the coii 
drivers can vary the currents. The size of sensing 
volume 50 within which the sensor will receive fields 
within the preselected field magnitude range from all 
coils will also depend upon the positions of the coils. 
However, for a typical system having three coils spaced 
at vertices of an equilateral triangle with sides about 
40 cm long, the sensing volume includes a region 
extending upwardly about 60 cm from the plane of the 
coils. At the plane of the coils, the sensing volume 



WO 96/41119 

PCT/US96/08411 

-16- 

r ils CT ■ 

5 location of o nl y one object «!L Position and 
ana plural sensors may be ut 1U «7 "T **«*" 

fiducial marher My include a s J 6 bo ^ 1 ;^" *«h 
10 and a sensor 54 of the SMB ^ ' tag 53 

connected to «T ^ " SenS ° r 18 "leasably 

connected to the sensor body. The amplification *Z 
converson unit 30 is ™„„ ^ ^ «P-mication and 

ju is connected to ^ 
additional sensors in the same manner „t 
The ta, of each sensor may be a ffixad to £ 

" Patlent " — « of the patient may bTacguir^ 

using any imaging method which win „ acqulrea 

desired imaging method. Where v '^ 91 " 9 ln **" 

the tag may be radioopagueT where Z 

the tag may include a ma^ia! Z,T^ ^ 
resonance response readily S^LTt? M9MtiC 

» bodily tissue. Ordinarily ^! T ° f 
. riJ - y ' the image is acquired 

to use of the catheter and probe discussed ™ 
sensors 54 typically are not present !„ «T 
step. Mter the imaging steo «1 1 * a9iB9 
sensors of the fiducial ^ «uT" T" " "* 

0 patienfs body at the positC of ^ ^n T 
catheter is employed, sensor siguals from sensors ™ 
the probe in the catheter are acoui-.,, 

sysT dT ^ ~ " - ~=-™r? ^ 
system determines the location and orf , 
5 54 k— . orientation of sensor 

54, and hence of fiducial markers 57 ,. • 

fields from coils 10 in the samT. ' 
determines the position anT or ^ " " ~ 
P"be 14. Xh. acguired position^ort^ £ £ 



WO 96/41119 



PCT/US96/08411 



-17- 

fiducial markers 52 can be used to register the 
previously acquired image of the patient with the 
position and orientation data for sensor 18 and probe 
14. For example, where the position and orientation 

5 data for the sensor 18 and probe 14 are displayed as a 
pictorial, image of the catheter tip 16 on a screen, the 
position and orientation of fiducial markers, as 
determined by the magnetic field locating method can 
also be displayed. The previously acquired image data 

0 may be displayed as a pictorial representation of the 
patient 1 s structures and that representation will also 
include a picture of the tags 53. The previously 
acquired images of the tags are brought into 
registration with the images of the fiducial markers 

5 derived from magnetic location by transforming one or 
the other image until the picture of the tags in the 
previously acquired image overlies the representation of 
the fiducial marker in the magnetically acquired image. 
Such registration can be achieved by manually adjusting 

0 inputs to an image display system while visually 
matching the displayed images, as described in the 
aforementioned "562 International Publication, or else 
by automatic calculation of the transformation and 
rotation parameters necessary to bring the images into 

5 registration. Such parameters can be calculated by 
equating the positions of three or more fiducial marker 
sensors in the magnetic location frame of reference (the 
X,Y,Z coordinates of Fig. 1) with the positions of the 
tags 53 in the image frame of reference as modified by 

» a matrix of unknown rotations a further matrix of 
unknown transformations. As further described in 
Appendix B, the resulting matrix equation yields an 
overdetermined system of nonlinear equations which 
include the rotation angles and transformation 

5 distances. When the images are properly registered, the 
position of the catheter tip is displayed in true 
location and orientation relative to the patient f s 
internal structures. In similar fashion, a greater 
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used on a given cycle may be adjusted to provide fields 
within the desired range of magnitude for one sensor. 
Signals from other sensor can be acquired during this 
cycle, although such other sensor may encounter a field 
5 outside the preselected range. In some cases, the coil 
currents are fixed for some cycles and adjusted for 
others. For example, where numerous fiducial markers 
are employed, it may be possible to achieve acceptable 
registration accuracy using fixed coil currents to 
10 provide the fields when data is acquired from the 
fiducial markers, whereas the coil currents may be 
adjusted in the manner discussed above when data is 
acquired from the sensor in the active device or 
catheter . 

15 Certain sensors tend to lose accuracy when 

exposed to magnetic fields above a predetermined 
maximum. For example, certain magneto-resistive sensors 
temporarily lose accuracy if exposed to magnetic fields 
above about 4 Gauss. Where such sensors are employed in 

20 a multi-sensor system, the system should have 
appropriate provisions to avoid exposing sensors to 
excessive fields. For example, if sensor 54a is 
disposed close to coil 10b, whereas sensor 18 is 
disposed remote from coil 10b, sensor 54a may be exposed 

25 to extremely high field levels if the system adjusts the 
current in coil 10b to the level required to generate 
the appropriate field at sensor 18. To avoid this 
condition, the data from all sensors can be acquired 
during all of the cycles, and control unit 34 can be 

30 arranged to increase coil currents in a progressive 
manner over several cycles when an increase is required. 
The data from sensors which are not used for position 
monitoring on a particular cycle can be used to inhibit 
further increases if the field at the nominally unused 

35 sensor is approaching dangerous levels. Thus, if the 
system is in the process of progressively increasing the 
coil current in coil 10b to provide an adequate field 
level at sensor 18, the system may terminate such 
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multiple carrier frequencies, equal to the number of 
sensors, the carrier frequencies for each coil being 
different from one another and different from all of the 
other carrier frequencies for all of the other coils. 
The components in the signals from the first sensor 
varying at the respective first carrier frequencies can 
be detected and the strength of the coil currents at the 
respective first carrier frequencies can be adjusted 
accordingly, so that the field magnitudes at the first 
carrier frequencies prevailing at sensor 18 will be 
within the desired range. In this stage of the process, 
the components varying at the other carrier frequencies 
are ignored. The process is reversed with respect the 
fiducial marker sensor 54a. This approach can be 
extended to greater numbers of sensors. Also, the 
frequency multiplexing and time multiplexing approaches 
can be combined. Thus, each coil can be driven at only 
one carrier frequency. During a first cycle, first 
sensor 18 is actuated and the components varying at all 
of the carrier frequencies are monitored. The coil 
currents are adjusted accordingly to bring the fields at 
each carrier frequency into the desired range. During a 
second cycle, the second sensor 54a (the first fiducial 
marker sensor) is used and the coil currents are 
adjusted in the same manner, and so on with the other 
fiducial marker sensors. 

In a further variant, the coil currents, and 
hence the strengths of the fields from the individual 
transmitters, can be adjusted to bring the field 
components detected by each individual component sensor 
to within a preselected range. In such a system, each 
coil is adjusted separately with respect to each 
component sensor. In a first cycle, the current to 
transmit coil 10a is adjusted to bring the individual 
sensor signal from component sensor 20, representing the 
field component magnitude in the local X' direction, 
into a preselected range. During this cycle, the 
signals from the other component sensors 22 and 24 are 
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sensor. The feedback control system thus operates to 
adjust the coil currents until the field component 
magnitude is at such single value. The position and 
orientation are calculated in the same way as discussed 

5 above. This variant has the advantage that nonlinearity 
in the component sensor response cannot affect the 
accuracy of the system. Provided that it is known that 
a particular reading from the sensor corresponds to the 
preselected value of field component magnitude, 

10 deviation from a linear relationship between component 
magnitude and sensor reading at other values of the 
component magnitude will not impair the accuracy of the 
system. 

Certain field sensors exhibit so-called "off- 

15 axis sensitivity". That is, the transfer function or 
relationship between field component magnitude along the 
sensitive axis of a particular component sensor and the 
reading from that component sensor varies when a strong 
field component orthogonal to such axis is present. 

10 Off-axis sensitivity can be corrected by using the 
readings from two component sensors to evaluate the 
magnitude of the field perpendicular to the sensitive 
axis of the third component sensor, and using that 
magnitude to determine a correction factor to be applied 

15 to the reading from the third component sensor. 

In the system described above, when the coils 
are driven with direct current, amplification and 
conversion unit 30 samples the data from each component 
sensor at a preselected delay time after the inception 

0 of the current. In a variant of this system, the 
command unit adjusts the delay time from inception of 
the current to sampling depending upon the current 
applied in a particular actuation. After inception of 
the current, the current gradually rises to its steady 

5 state value, and hence the overall field strength or 
dipole moment likewise rises gradually. The changing 
field induces eddy currents in electrically conductive 
materials which may be present in or near the sensing 
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example, the transmitting coils 10 of Fig. i can be 
disposed in a vertical plane arrangement as illustrated 
for the sensors in Fig. 5. Also, the number of 
transmitters and sensors can be varied. For example, as 

5 disclosed in the aforementioned International 
Publication WO94/04938, the system may include a sensor 
having three mutually orthogonal sensing axes and a set 
of three mutually orthogonal transmitting coils having a 
common center in the fixed frame of reference. 

0 Alternatively, the system may include a single uniaxial 
transmitting coil on the object to be located and three 
sets of sensors on the affixed frame of reference, each 
such sensor set incorporating three receiving coils or 
component sensors adapted to detect fields in a mutually 

5 orthogonal direction. The reciprocal arrangement — 
with three sets of mutually orthogonal transmitting 
coils and a single uniaxial sensor on the object to be 
tracked — is also usable. In general, the sensors and 
transmitters should define a plurality of transmitter- 

0 receiver pairs, each including one element on the object 
and one element on the fixed frame of reference. 

In a further variant, the system adjusts the 
output of the transmitters in response to the calculated 
disposition of the object being tracked, rather than 

5 directly in response to the component signals or total 
field magnitude signal. Thus, the system can initially 
operate with the default of current values; derive an 
initial reading of the position and orientation of the 
object and then use that initially determined position 

3 and orientation to calculate the desired setting for 
each coil to achieve the desired field levels at the 
sensor. Such desired setting is selected to yield the 
field within the desired magnitude range at the object, 
assuming that the object has the position and 

5 orientation found in the, initial reading. On the next 
cycle, the so-calculated coil currents are utilized and 
the process is repeated. In a variant of this approach, 
the system can store a lookup table listing the 
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appropriate coil currents for 
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Appendix A: calculation of Position and Orientation 

5 Provided that we know the physical 

configuration of the field generator (transmitter) 
which is fixed in position during operation, the 
magnetic field detected by a sensor is a function of the 
position and orientation of the sensor. in our system 
10 field generator coils are stimulated sequentially. The 
field sensed by a sensor (3 component sensors per probe) 
can be expressed in terms of position in x, y, z and 
orientation a, P, y (roll, pitch and yaw respectively) 
i.e.: 

15 

B[ sensor] [coil] = f [sensor] [coil] (x, y, 

z, a, p, y) 

Where [sensor] designates a particular sensor 
and [coil] designates a particular transmitter coil. 
20 If the real field the sensor measures when 

[coil] is on is B« [sensor] [coil] , then theoretically, 

B» [sensor] [coil] = B[ sensor] [coil] 

25 

i.e. 

B» [sensor] [coil] - f [sensor] [coil] (x, y, 
z, a, p, y) = o.o 

30 

Since we have 3 sensors and 3 coils, the total 
equations are 9 with 6 unknowns (x, y, z for probe space 
location, a, P, y for its orientation). By applying non- 
linear least square method, we can find unique x, y, z, 
35 a, P, y for the probe. 

Above shows the general ideal of the 
algorithm. In detail: 
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is: 



10 



15 



em e P n e,\i 
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e P 3i e,n e,3 3 



and, since the threp can. 

be orthogonal to each ot^r their ** ^ *** "* 
can described as: * non - ort »°9°nal axes 
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10 



since we use roll (a), pitch (P) , yaw(y) to 
define probe orientation, ortho_OV[i] [ j] can be also 
described by: 



sin(P) sin(y) 
sin(P) cos(y) 



15 sin(P) sin(y) 
sin(P) cos(y) 



ortho_OV[l] [1] = cos(a) cos(y) - sin(a) 

ortho_OV[l][2] o cos(a) sin(y) - sin(a) 

ortho_OV[l][3] = -sin(a) cos(p) 

ortho_OV[2][l] = -cos(P) sin(y) 

ortho_OV[2] [2] = cos(P) cos(y) 

ortho_OV[2][3] - sin(P) 
ortho_OV[3][l] - sin(a) cos(y) + cos(a) 

ortho_OV[3)[2] = sin (a) sin(y) - cos (a) 

ortho_OV[3] [3] - cos(a) cos(p) 



20 orthogonal vector inatrix can therefore be 

calculated by matrix multiplication of previous defined 
matrix T and orhto_OV: 

ov=T*ortho_OV 

25 

The theoretical magnetic field for an 
orthogonal system generated at sensor position pointing 

to direction can be expressed as: 

30 f[coiT\[i\(x >y ,z,a,p t y) 

(detail abbreviated), where / is a known function, and 
includes a dipole moment term having magnitude 
proportional to the current flow in the particular coil. 
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»ay not J^Z^T" 9 "* 11 ** °°™ e "<» (— «. 
y nor toe Perpendxcular to each oth«*.\ 

field sensor measures shouid be ' M,netiC 



Assu,. the real field sensor detected when 
10 coil is on is B lsensor][coil] 

i tnen: 

i-tannx-HoX/] - B(samr][colt\ « 0.0 

therefore, the 9 equations to be solved for x 
15 y. z, a, p, r , are. 
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B'[sensor2][coil2]- Jl/[coil2m(x,y,z,a,fi t r) x ov -\i][sensar2] 



M =0.0 



3 

B'[sensor2][coil3}- ^f[coil3][i](x,y,z,a,fi,r) * ov' l V][semor2) = 0.0 

5 



B'[semor3)[com]-Jlf[com][i](x,y,z,a,J3,r) x ov '\i][sensor3] = 0.0 



B , [sensor3)[coil2]-Y l f[coil2][i](x,y,z,a,fi,r) x ov"'[0[«n«w3] = 0.0 



i=i 



10 

3 

B , [sensor3][com]- ^£f[coil3][i](x,y,z,a,fi,r) * ov'^*™^] = 0.0 



i=i 



A well-known non-linear least squares equation 
15 solver can be applied to solve the above equations and 
find probe position x, y, z and orientation a, p, y. 
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Appendix B: Registration of Fiducial Markers 

Assumption: 

WS knoW the co °rdinates of three fldn^.i 

» Hy3. ((MX1 ' Myl ' "">' "«>. (—.' 



IS 



and the coordinates of the ^ 
fiducal points in the X Y Z ™« „ \ ***** 
fixed , ' coordlnat e system of the 

nxed frame of reference of the 

apparatus: magnetic locating 



Py3, Pz3)} 



{(PX1, pyx, P21) , (px2/ ^ pz2)f ^ 



20 



25 



-» ------ - 

Tne rotation matrix R is: 



i?nc iter 
toy Ryy Rzy 

_Rxz Ryz Rzz 



30 



where 
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Rxx = cos(a) x cos(y) - an(a) x an(fi) x an(y) 
Rxy = cos(a) x anOO + sin(a) x x coh» 
Rxz = - sn(or) x cos(fi) 
Ryx = -cos(J}) x an(y) 
J?>y = cosO?) x cosOO 
Ryz = an(/9) 

Jfcx = an(a) x coa» + cos(a) x an(£) x anOO 
Rzy = an(a) x anOO - cos(a) x an(0) x coa» 
Rzz = cos(a) x cos(^) 



The translation matrix T is: 



X ^ Z 

x y z 
x y z 



The relationship between image position and 
magnetic location apparatus frame of reference position 
is: 



10 



or 



'Mxl 


Myl 


Mzl' 




~Pxl 




Pzl 








Mx2 


Myl 


Mzl 


R+T = 


Pxl 


Pyl 


Pzl 








Mx3 


Myl 


Mz3_ 




Px3 


PV3 


Pz3 








'Mxl 


Myl 


Mzl 




'Pxl 


Pyl 


Pzl 




"0 


0 


0" 


Mxl 


Myl 


Mzl 


R + T- 


Pxl 


Pyl 


Pzl 




0 


0 


0 


Mx3 


My3 


Mz3_ 




Px3 


Py3 


Pz3 




0 


0 


0 
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This forms 9 single equations: 
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Mx 1 x Jbac + My 1 x Rxy + Mz 1 x Rxz + * - Px \ = 0 

Airlxi2zr + M>-lxi^ + ^ lxite + z _ Pzl = 0 
x Rxx + Kfy2 x io? + A*2 x /6s + x - Px2 = 0 

Mx2 xtex + My2xR Z y + M22xRz2+z -p 2 2 = 0 

Mx ^Ryx + My3xRyy + M 23xByz+y _ Py2 = 0 
M x3xtex + My3xR Z y + Mz 3xR zz + 2 - Pz3 = 0 



(*) 



with 6 unknowns: a, B# v for r«f S f< 
' V, 3 for sp.ce tr„,l.tion. ^ M9laS ' 

" 9iStratiCn *~ *« *• ^ rotation 
the above 9 equations. "xving 

L0 

in pseudo-code, the registration procedure is: 
registration () 

Initialize data buffer* 

initialize magnetic locating system; 

0 fiducial points;™ ^ °* «" — 

l ea «+ APPly equation sol ving method (non-linear 

least square) on equations(*) , * r 

and find rotation angles a, B v and 
space translation x, y, 2 ; and 

POrm rotation ^trix R from a, fi T and 
translation matrix T from x f y, 2 ; 
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The rotation matrix R and translation matrix T 
can be applied to subsequent positions obtained, 
including the positions of the of the sensor associated 
with the catheter. 



WO 96/41119 



-36 



PCT/US96/D8411 



10 



IS 



20 



25 



CLAIMS: 

1. Apparatus for determining position 
comprising: 

(a) field generating means ( i 0 , for 
Producing a plurality of M g netlc £i . laa _ ^ ^ 
fxeld having parameters varying with location within a 
sensing volume (5 o, according to a known pattern of 
variation, said patterns of variation heing different 
for dxfferent ones of said fields; 

(b) at least one sensor (i 8 ) adapted to 
detect one or .ore of said parameters prevailing at such 
sensor when said sensor is at an unknown location within 

said sensing volume and to provide one ™- 

F viae one or more sensor 
signals representing said one or more parameters; 

(=) calculation means (32 ) for 
aetermining a calculated location of said sensor based 
upon said sensor signals and said Known patterns of 
variation; and 

(d) feedback control Means (34) for 
adjusting said field generating means to alter said 

known pattern of variation of at 

1 at -Least one of said 

fields responsive to said 

said sensor signals, sa id 

calculated location or hoth so as to maintain said 
parameters of each such altered fi.i a prevailing at said 
sensor within a preselected range. 

2. Apparatus for determining position 
comprising: 
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(a) a plurality of sensors (110) adapted 
to detect one or more parameters of a magnetic field 
prevailing at each such sensor and to provide one or 
more sensor signals representing said one or more 

5 parameters, said sensors being disposed adjacent a 
sensing volume (150) at different positions, different 
orientations or both; 

(b) field generating means including at 
least one magnetic field transmitter (118) movable to an 

10 unknown location within said sensing volume and adapted 
to provide at least one magnetic field in said sensing 
volume, so that the magnetic field from each transmitter 
varies in according to a known pattern of variation with 
location relative to such transmitter whereby said 

15 sensor signals will represent one or more parameters of 
the field or fields from said at least one transmitter; 

(c) calculation means (32) for 
determining a calculated location of said at least one 
transmitter based upon said sensor signals and said 

20 known patterns of variation; and 

(d) feedback control means (34) for 
adjusting said field generating means to alter said 
known pattern of variation of at least one of said 
fields responsive to said sensor signals, said 

25 calculated location or both so as to maintain said 
parameters of each such altered field prevailing at each 
said sensor within a preselected range. 
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3. 



Apparatus as claimed in claim ! or claim 
2 wherein said field generating means is operative to 
Provide DC magnetic fields within said sensing volume. 

4. App aratus as claiMd in oiain 3 Kherein 

sa,d field generating means includes a pluraUty Qf 
coiis (10 . i 2 „, 122 , 124) anfl Mans intermi 
actuating each said coil (34 , 42 , 44 , 46) wlth , 
current, said feedhac* control means inciudes means for 
altering the magnitude of each said direct current, and 
wherein said calculation means includes maans tor 
sampling the signals £ro m each sensor when a delay time 
has elapsed after commencement of each said intermittent 
actuation, said f.edoacK control means further including 
means for varying the del ay «« for 

accordance with the magnitude of the current applied in 
such actuation so as to reduce the delay time when 

lesser currents are employed ar.a s — 

employed and increase the delsy time 

"hen greater currents are employed. 

5. Apparatus as claimed in claim i or claim 
* wherein said field generating means is operative to 
provide AC magnetic fields within said sensing volume. 

«• Apparatus as claimed in claim 2 or elaln 
wherein each said sensor has merimum accuracy within a 
preselected optimum operational range for each said 
parameter and said preselected range for each said 
Parameter suhstantially corresponds to said optimum 
operational range. 

7. Apparatus as claimed in claim i or claim 

2 wherein said preselected range for each ^ 

ye lor each said parameter 
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consists of a single value, and wherein said feedback 
control means is operative to maintain each said 
parameter of the field at each sensor substantially at 
said single value. 
5 8. Apparatus as claimed in claim 1 wherein 

said field generating means includes a plurality of 
transmitters (10) disposed adjacent said sensing volume. 

9. Apparatus as claimed in claim 2 or claim 
8 wherein each said transmitter is operative to radiate 

10 a field having magnitude decreasing with r n where r is 
distance from the transmitter and n is a number greater 
than l. 

10 . Apparatus as claimed in claim 9 wherein 
each said transmitter is operative to radiate a field 

15 having magnitude substantially equal to K/r 3 where K is 
a real number, and wherein said feedback control means 
is operative to adjust K. 

11. Apparatus as claimed in claim l further 
comprising a sensor body (18), said at least one sensor 

10 including a plurality of component sensors (20, 22, 24) 
disposed on said sensor body, each said component sensor 
being operative to measure the magnitude of a magnetic 
field component in a preselected local direction 
relative to the sensor body. 

15 12. Apparatus as claimed in claim 11 wherein 

said field generating means includes a plurality of 
transmitters (10) disposed at spaced-apart locations 
adjacent said sensing volume, said feedback control 
means being operative to adjust each said transmitter to 
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trans.it a plurality of fleloE eo ^ ^ ^ ^ 
field, the magnitude of a component of such field i„ 
the local direction of one said component sensor will be 
within a preselected range. 

13. Apparatus as claimed in claim n wherein 
each said component sensor is a Hall effect sensor, . 
magnetoresistive sensor, a magnetooptical sensor or a 
flux gate magnetometer. 

14. Apparatus as claimed in claim 1 or clai, 
2 wherein each said sensor includes a ceil and wherein 
said field generating meens is arranged to provide AC 
fields in seia sensing volume. 

15. Apparatus as claimed i„ claim x Kherein 
said field generating means includes . p^.^y of 

» transmitters disposed at sp.ced-.pert lections en one 
side of . reference plane and wherein seid sensing 
vciume extends to the other side of seid reference 
plane. 

16. Apparatus es claimed in cleim 15 wherein 
seid transmitters ere disposed substantia!^ i„ a common 
plane. 

17. Apparatus as claimed in claim 2 Kheretn 
seid pluraUty of sensors are disposed et speed-apart 
lecations on one side of a reference plena and wherein 
said sensing volume extends to the ether side of said 
reference plane. 

«• Apparatus as claimed in cleim i 7 
said sensers are disposed substantially m a c^men 
plane. 
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19. Apparatus as claimed in claim 1 or claim 
2 wherein said at least one sensor includes a plurality 
of sensors, and wherein said sensors can be impaired by 
exposure to fields having components above a 
5 predetermined threshold magnitude, and wherein said 
feedback control means includes means for controlling 
said field generating means to assure that no sensor is 
exposed to a field component having magnitude above its 
threshold value. 

10 20. Apparatus as claimed in claim 1 or claim 

2 wherein said at least one sensor includes a plurality 
of sensors and wherein said sensors can be impaired by 
exposure to fields having components above a 
predetermined threshold magnitude but can be restored 

15 after such exposure upon application of a reset routine, 
said feedback control means including means for 
initiating such reset routine in the event of such 
exposure . 

21. Apparatus for sensing the disposition of 
20 an object in a frame of reference comprising: 

(a) transmitter means (10) including at 
least one transmitter for providing a field; 

(b) sensor means (18) including at least 
one sensor for sensing one or more parameters of a field 

25 prevailing at each such sensor and providing one or more 
sensor signals indicative of such parameters, said 
transmitter means and said sensor means cooperatively 
defining a plurality of transmitter-sensor pairs, each 
including one transmitter and one sensor as elements of 
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the p. ir , C na element of Mch „ uch pair ^ 
on the object and the other element of each such pair 
oein, disposed at a known disposition in said frame of 
reference ; 

(c) calculation mea ns (3 2) for 
determining the disposition of said object in said fratte 
of reference based upon said sensor signals; and 

(d) feedback control means (34) for 
adjusting said transmitting means responsive to said 

" sensor signals to thereby maintain at least one of said 
sensor signals within a preselected range. 

22. a method of determining the disposition 
of an object in a frame of reference comprising the 
steps of: 

5 Providing one ot- 

* wne or more magnetic 

"elds eaoh having a known pattern of variation with 

respect to said frame of reference; 

(h) sensing one or more parameters of 

saxd one or more magnetic , ields prevailing at said 
- object to provide one or more sensor signals indicative 

of said one or more parameters; 

(o) calculating the disposition of said 

object in said frame of reference hased upon said sensor 

signals and said known pattern or patterns of variation- 
25 and ' 

(d) altering said known pattern of 

variation of at least one of said one «• ™« 

aia one or more magnetic 
fields responsive to said sensor signals, said 
calculated disposition or both so as to maintain said 
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parameters of each such altered field prevailing at said 
object within a preselected range. 

23. A method of determining the disposition 
of an object in a frame of reference comprising the 
5 steps of: 

(a) providing one or more magnetic 
fields each having a known pattern of variation in space 
with respect to object; 

(b) sensing one or more parameters of 
to said one or more magnetic fields prevailing at one or 

more sensors having known disposition in said frame of 
reference to provide one or more sensor signals 
indicative of said one or more parameters; 

(c) calculating the disposition of said 
15 object in said frame of reference based upon said sensor 

signals and said known pattern or patterns of variation 
of said one or more magnetic fields; and 

(d) altering said known pattern of 
variation of at least one of said . one or more magnetic 

10 fields responsive to said sensor signals, said 
calculated disposition or both so as to maintain said 
parameters of each such altered field prevailing at at 
least one of said sensors within a preselected range. 
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